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1. INTRODUCTION 

The lack of a consistent quantum theory of the full gravitational field and its sources has 
stimulated the development of quantum cosmology, a less complete but more tractable method 
to investigate the influence of quantum effects on the evolution of the universe. The primordial 
universe, when presumably curvatures and densities approach the Planck scale, is believed to be the 
privileged scenario in which the quantum aspects of gravity are expected to become important or 
even dominant. In its broadest sense, quantum cosmology consists in "freezing out" all but a finite 
number of degrees of freedom of the the gravitational field plus its sources (through imposition 
of symmetry requirements) and then quantizing the remaining ones. This procedure, initiated by 
DeWitt [1], is expected to provide some general insights on what an acceptable quantum theory of 
gravity should be like, although it cannot be strictly valid and is open to criticism [2]. Such a line 
of attack has been extensively explored to quantize Friedmann-Robertson-Walker(FRW) universes 
with varying matter content such as a scalar field [3-6], radiation [7-9], a spinor field [10], dust 
[9,11-14 ] or a Rarita-Schwinger field [15]. 

The present paper is dedicated to a further study of the quantum theory of a radiation-filled 
FRW universe. Differently from previous investigations of this system [7-9], here we discuss both 
closed and open models, deal only with normalizable wave functions, and pay full attention 
to the domain of self-adjointness of the Hamiltonian operator. We follow the Arnowitt-Deser- 
Misner(ADM) genuine canonical quantization method [16]. In this approach one has to solve the 
constraint equations at the classical level and go over to a reduced phase space spanned by inde- 
pendent canonical variables alone, and this process demands a definite choice of time. Although 
often leading to complicated and time-dependent Hamiltonians, this formalism has the great ad- 
vantage of reducing the problem to one of standard quantum mechanics, enabling one to make full 
use of the powerful theory of linear operators in Hilbert space. In our treatment the time variable 
is chosen as conformal time, as this enormously simplifies the form of the Hamiltonian operator, 
making the quantum dynamics exactly soluble. Further reasons for choosing conformal time are 
given in [8]. 

After the ADM reduction of phase space, only one degree of freedom remains, which is taken 
to be the scale factor R. Since R is restricted to positive values, it becomes necessary to impose 
boundary conditions on the wave functions belonging to the domain of the Hamiltonian operator 
to ensure its self-adjointness. For the simplest of such boundary conditions the time evolution of 
wave packets of Gaussian type is worked out, and it is shown in full generality for the first time 
that both the closed and open models are nonsingular. It is remarked that in the context of the 
ADM quantization the so-called tunneling boundary condition on the universal wave function is 
in conflict with self-adjointness of the Hamiltonian operator, or, equivalently, with unitarity of the 
quantum evolution. This is not an artifact of the particular quantization scheme adopted here, 
since the ADM and Wheeler-DeWitt descriptions are equivalent for the model at hand [7]. An 
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initial condition such that the probability density is sharply concentrated at R = and under 
which the classical regime sets in is considered. In the closed case, under such an extreme initial 
condition the probability of finding any finite radius for the universe vanishes. Therefore a density 
parameter Q = 1 becomes a prediction of the model, without neither sacrificing the requirement 
of square-integrability on the wave functions nor imposing boundary conditions inconsistent with 
the space of admissible state vectors . A physically questionable aspect of the initial condition 
adopted is pointed out, however. Stationary quantum geometries are shown to exist for the closed 
model, but the existence of oscillating coherent states of the geometry is precluded by the boundary 
conditions required to enforce self-adjointness of the Hamiltonian operator, a result at variance 
with previous findings [9]. 

The layout of this paper is as follows. In Section 2 the classical model is specified and the ADM 
reduction of phase space [7] is briefly reviewed. In Section 3 the problem of the necessary boundary 
conditions to ensure self-adjointness of the Hamiltonian operator is considered, and the respective 
propagators are written down when the two simplest of such boundary conditions are adopted. In 
Section 4 the motion of wave packets obeying different boundary conditions is obtained in closed 
form, and it is verified that the singularity is avoided in all cases. A special initial condition under 
which the quantum model is forced into the classical regime is discussed, and the consequences for 
the case of a closed universe are considered. Stationary quantum geometries are taken up in Section 
5, with particular emphasis on how the restricted domain of self-adjointness of the Hamiltonian 
operator influences such states. Section 6 is devoted to some final comments. 

2. DYNAMICS OF THE CLASSICAL MODEL 

The line element for a homogeneous and isotropic universe can be written in the FRW form 

(we take c = 1) 



where denotes the metric for a 3-space of constant curvature k = +1,0 or —1, corresponding 
to spherical, flat or hyperbolic spacelike sections, respectively. 

The matter content will be taken to be a perfect fluid, and Schutz's canonical formulation of 
the dynamics of a relativistic fluid in interaction with the gravitational field will be employed [17]. 
The degrees of freedom ascribed to the fluid are five scalar potentials tp, a, (3, 6, S in terms of which 
the four-velocity of the fluid is written as 



ds 2 = g vX dx v dx x = -N(t) 2 dt 2 + R(t) 2 a ij dx i dx j , 



(2.1) 




(2.2) 
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where \i is the specific enthalpy. By means of the normalization condition 

g vX U v U x = -l (2.3) 

one can express /i in terms of the velocity potentials. The action for the gravitational field plus 
perfect fluid is 

S=[ d 4 xy/=g~ ^R + 2 [ d 3 xVhh l7 K ij + / d 4 x^g~p (2.4) 

JM JdM JM 

in units such that c = 16nG = 1. In Eq.(2.4) p is the pressure of the fluid, is the scalar 
curvature derived from the spacetime metric g u \, hij is the 3-metric on the boundary DM of the 
4-manifold M, and K lj is the second fundamental form of the boundary [18]. The surface term is 
necessary in the path-integral formulation of quantum gravity in order to rid the Einstein-Hilbert 
Lagrangian of second-order derivatives. Variations of the pressure are computed from the first law 
of thermodynamics. 

Compatibility with the homogeneous spacetime metric is guaranteed by taking all of the velocity 
potentials of the fluid as functions of t only. We shall take p = (7 — l)p as equation of state for the 
fluid, where 7 is a constant and p is the fluid's energy density (we shall eventually put 7 = 4/3). 
In the geometry characterized by (2.1) the appropriate boundary condition for the action principle 
is to fix the initial and final hypersurfaces of constant time. The second fundamental form of the 
boundary becomes = —hij/2N. As described in its full details in [7], after inserting the metric 
(2.1) into the action (2.4), using the equation of state, computing the canonical momenta and 
employing the constraint equations to eliminate the pair (0,p e ), what remains is a reduced action 
in the Hamiltonian form 

S r = J dt{Rp R + <p P(p + S Ps - NH} (2.5) 

where an overall factor of the spatial integral of (deta) 1 ^ 2 has been discarded, since it has no effect 
on the equations of motion. The super-Hamiltonian Ti is given by 

H = - (j^ + 6kRj + pi R- 3 ^V e s . (2.6) 

The lapse N plays the role of a Lagrange multiplier, and upon its variation it is found that 
the super-Hamiltonian 7i vanishes. This is a constraint, revealing that the phase-space contains 
redundant canonical variables. 
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According to the ADM prescription, in order to perform a bona fide canonical quantization one 
must go over to a reduced phase space spanned by independent canonical variables alone. This can 
be achieved by first making a choice of time and then solving the super-Hamiltonian constraint 
equation TC = for the canonical variable conjugate to the time chosen in the first step. This 
ensures that the final action preserves its canonical form, and the Hamiltonian in the reduced 
phase space is identical to the canonical variable whose Poisson bracket is unity with whatever was 
chosen as time, but now expressed as a function of the remaining independent canonical variables 
[16]. In the conformal-time gauge N = R (that is, henceforward t denotes conformal time) and for 
7 = 4/3 (radiation) this procedure leads to the very simple reduced action [7] 

S r = J dt\kp R - (j| + QkR 2 ^ I . (2.7) 

Only one degree of freedom is left, namely the scale factor R, and the Hamiltonian in the reduced 
phase space is 

2 

H = ^- + 6kR 2 . (2.8) 
24 

Hamilton's equations of motion lead immediately to 

R + kR = . (2.9) 

The solution for R(t) can be written as 

!sint if k — +1 
t if k = (2.10) 
sinht if k — — 1 

with a suitable choice for the origin of conformal time t. The standard cosmic time r is related to 
conformal time by 

dr = Rdt , (2.11) 



hence 



1 — cos t if k — 1 
t 2 /2 if k = 
cosht — 1 if k — — 1 



(2.12) 



with the convention that r = when t — 0. In the spatially flat case (k — 0), for instance, one 
recovers the usual behavior R = C t 1 / 2 for the scale factor [19]. It is seen that Hamilton's principle 
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based on the reduced action (2.7) gives rise to the same equations of motion as those obtained 
by first varying the full action (2.4) and then simplifying them through the use of the spacetime 
symmetries of homogeneity and isotropy and of the equation of state p = p/3 . Such a consistency 
check is indispensable if quantization in minisuperspace is to have any meaning at all. 



3. QUANTIZATION, SELF-ADJOINTNESS AND BOUNDARY CONDITIONS 

The remarkably simple form of the Hamiltonian (2.8) makes it possible to find exact results for 
the cosmic evolution at the quantum level. The quantum dynamics is not so straightforward as 
one might think at first sight because the scale factor R is restricted to the domain R > 0, so that 
the minisuperspace quantization in the i?— representation deals only with wave-functions defined 
on the half-line (0, oo). It is well-known that in such circumstances one usually has to impose 
boundary conditions on the allowed wave functions otherwise the relevant operators will not be 
self-adjoint, the most important of all operators being the Hamiltonian, that must be self-adjoint 
in order that the time evolution be unitary. The need to impose boundary conditions to ensure self- 
adjointness has been long recognized by practitioners of the ADM formalism as applied to quantum 
cosmology [3,11]. Very recently it has also been seen to have non-trivial cosmological implications 
in the Wheeler-DeWitt approach [20]. What does not appear to have been duly emphasized is 
that self-adjointness conditions depend on the set of allowed state vectors. It has been argued by 
adherents of the many-worlds interpretation [8] and of the pilot-wave formulation [21] of quantum 
theory that non-normalizable wave functions are unavoidable in quantum cosmology. Irrespective 
of whether their arguments are physically well founded or not, what we want to stress here is that 
the self-adjointness conditions must be consistent with the point of view adopted, if the results 
obtained are to be regarded as trustworthy. 

As follows from the substitution p R — > —id/dR, the Hamiltonian operator associated with the 
classical Hamiltonian function (2.8) is (we take h = 1) 

H -T^ + 6kR2 <"> 

defined on the half-line (0, oo). The condition for H to be symmetric (which, in turn, is a necessary 
condition for H to be self-adjoint) is 

^M 2 ) = {H^M (3.2) 



or 
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poo rl^ihr, r°° rPih* 



where the asterisk stands for complex conjugate. Integrating by parts twice this leads to 



If both ip and its derivative are square-integrable, the left-hand side of (3.4) vanishes and we 
are left with 

(^-f, 2 )(0H0. (3,, 

Then it can be shown [22] that to ensure the validity of this condition it is necessary and sufficient 
that the domain of H be restricted to those wave functions such that 

^'(0) = aip(0) (3.6) 



with a G (— oo, oo]. This generic boundary condition was explicitly taken into account in [11] and 
implicitly used in the simplest cases a = and a = oo in [7], in which only square-integrable wave 
functions were considered acceptable. 

If the potential is unbounded from below the Hamiltonian may possess eigenstates ip G L 2 (0, oo) 
such that ip' is not square-integrable. In this case, or if non-normalizable wave functions are allowed, 
the correct condition for symmetry of the Hamiltonian is Eq.(3.4), and as such it has been recently 
employed in [20]. However, in [8] non-square-integrable wave functions were argued to be necessary 
in quantum cosmology but, inconsistently with this point of view, Eq.(3.6) was imposed on the 
allowed wave functions to allegedly enforce self-adjointness of the Hamiltonian operator. As a 
consequence of demanding that the initial wave function when t = be perfectly localized at the 
singularity R = 0, Tipler [8] was led to a universal wave function of the form 



iP(R,t) 



3i 



ALp sint 



1/2 

2 



cxp 



(3tt /4i) (cot t){R/Lpf = A{t)exp[iB{t)R 2 } , (3.7) 



where L P denotes the Planck length, and which satisfies the boundary condition (3.6) for a = 0. 
Since this wave function is not square-integrable, it should obey Eq.(3.4), that is 

lim UiA*BR] = , (3.8) 

R^oo 



8 



which is not satisfied because both A(t) and B(t) are different from zero. Therefore, a wave function 
of the form (3.7) cannot represent a possible state of the universe. This means that Tipler's initial 
condition on the wave function of the universe together with his imposition of boundary condition 
(3.6) with a = are in conflict with the requirement that H be a self-adjoint operator, and this 
renders the conclusions of [8] invalid. 

In the present work we shall deal only with square integrable wave functions, so that the set of 
admissible states in the i?- representation is the Hilbert space L 2 (0, oo). Therefore, the domain of 
self-adjointness of the Hamiltonian operator is restricted to those wave functions that obey (3.6). 
For the sake of simplicity, here we shall address ourselves only to the cases a = and a = oo, that 
is, the boundary conditions we shall be concerned with are 

^'(0,f) = (3.9a) 

or 

ip(0,t) = . (3.96) 



Both of these conditions refer to what happens to a wave packet when it hits the singularity 
R = 0. The boundary condition (3.9b) was advocated by DeWitt to keep wave packets away from 
the singularity, but in general it is not powerful enough to prevent wave functions from becoming 
concentrated in the neighborhood of R = [23]. As a matter of fact, it has been argued [11] that 
DeWitt's boundary condition is just not relevant to the issue of quantum gravitational collapse . 

The time development of the models is fully determined once one is in possession of the prop- 
agator or Green's function. Let G(x,y,t) be the propagator for the problem in the usual Hilbert 
space L 2 (— oo, oo). Then the propagator for the problem in the restricted Hilbert space L 2 (0, oo) 
is 

G (a) (R, R', t) = G(R, R', t) + G(R, -R', t) (3.10) 

if the boundary condition is (3.9a), or 

G ib) {R,R',t) = G(R,R',t)-G(R,-R',t) (3.11) 

if the boundary condition is (3.9b), as noted by several authors [24,8,25]. 
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The general Green's function for the Hamiltonian (3.1) on the usual Hilbert space L 2 (— oo, oo) 

is 



G(x,y,t) 



in sin 



(Vkt) 



1/2 , „. rr 

expj s . n( ^ t) [(x 2 + V 2 ) cos(Vkt) - 2xy] j , (3.12) 



as one immediately obtains from the expression of the propagator for the harmonic oscillator [26] 
by setting m — 12 and uo = \fk. In the limiting case k = the well-known free-particle propagator 
is regained, whereas for k = — 1 all one has to do is make use of the simple formulae cos(it) = cosht 
and sin(ii) = isinht. The latter case corresponds to the quantum mechanics of a particle in an 
inverted oscillator potential, studied extensively in [27]. 

In the quantum cosmology a la Hartle-Hawking-Vilenkin-Linde an essential role is played by 
initial or boundary conditions. One of these is the so-called tunneling boundary condition [28,29], 
according to which the wave function of the universe must consist only of outgoing modes at 
singular boundaries of superspace. In our present context this amounts mathematically to 

where J is the probability current density. However, from (3.6) it follows immediately that J = 
because a is a real number. One is forced to conclude that, at least for this minisuperspace 
model, the tunneling boundary condition cannot be implemented because it is irreconcilable with 
self-adjointness of the Hamiltonian operator or, equivalently, unitarity of the time evolution. Fur- 
thermore, since the ADM and Wheeler-DeWitt descriptions are equivalent for the present model 
[7], this difficulty is not an artifact of the particular quantization scheme adopted here. 



4. EVOLUTION OF THE QUANTUM MODELS 

We shall now dedicate some paragraphs to the description of the main features of the dynamical 
evolution of our quantum cosmological models. This will be done by first following the time 
development of wave packets and then by studying the effect of imposing a very special initial 
condition. 

4.1 Motion of Wave Packets 

Let us start by working out the dynamical evolution of representative initial wave packets. The 
first initial state to be considered is the one described at t = by the normalized wave function 
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(4.1) 



where j3 = a + ip with p real and er > 0, corresponding to the boundary condition (3.9a). 
The initial wave function (4.1) is an even function of R, so that 

(R, t) = / G (a) (R, R', t)4 a) (# W 

JO 

P OO /* OO /* OO 

= / G{R,R',t)^ a) (R')dR! + / G(R,-R',t)^ a) {R')dR' = / G{R,R',t)^ a) (R')dR' . (4.2) 

JO JO J — oo 



Inserting the propagator (3.12) and the initial wave function (4.1) into (4.2) and performing the 
Gaussian integration one finds 

= (^) 1/4 (^= _} 1/2 

K n J \cos(Vkt)[(3tan(Vkt) - 6iVk] J 



x exp< 



tan( v / A; i) \ cos 2 (Vk t) \j3 taa(Vk t) - 6iVk] ) 



(4.3) 



An important quantity is the expectation value of the scale factor 



(R)l a) = f°° R\^«\R,t)\ 2 dR , 
Jo 



(4.4) 



which can be readily computed from the wave function (4.3). We find 



{R)i a) 



1 / 2 



<7 2 sin t + (6 — ptant) 2 cos 2 1 if k — +1 



a 2 t 2 + (6 - pt) 2 



if fc = 



(4.5) 



ya 2 sinh 2 t + (6 -ptanht) 2 cosh 2 t if k = -1 

Notice that, in all cases, (R)^ never vanishes. For A; = 0orA; = — 1 and p > the expectation 
value of the scale factor initially decreases, reaches a minimum value and then grows steadily 
without limit, whereas if p < there is a continuous expansion without bound. As expected, 
for large t the highest expansion rate belongs to the hyperbolic model (k = —1). For k = 1 the 
universe oscillates between a minimum and a maximum radius. An interesting interpretation of 
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this behavior in terms of reflection of parts of the wave packet as they hit the origin R = can be 
found in [11]. 

It should be stressed that the previous results establish that, at the quantum level, the singular- 
ity is avoided in all cases (that is, k — 0, ±1) according to the following reasonable criterion [11,30]: 
the quantum system is singular at a certain instant if (il)\f\il>) = for any quantum observable / 
whose classical counterpart / vanishes at the classical singularity, ip being any state of the system 
at the instant under consideration. For FRW models the relevant quantum observable is / = R, 
since R — defines the classical singularity. This criterion is in consonance with the usage in 
quantum cosmology. Indeed, since R is a positive operator on L 2 (0, oo), if (R) t = then ip(t) is 
sharply peaked at R = 0, and a strong peak in the wave function at a certain classical configuration 
is regarded in quantum cosmology as a prediction of the occurrence of such a configuration [29]. 

We now turn our attention to the boundary condition (3.9b). As initial wave function let us 
choose 



128a 3 



1/4 



Re~^ 2 



IT 



(4.13) 



Taking advantage of the odd character of this wave function, we can write 



P OO P oo 

(R, t) = / G {b) (R, R', t)4 b) (R')dR' = / G{R, R', t)^ (R')dR' 

JO J-oo 

Insertion of (4.13) and (3.12) into (4.14) yields 



128(7 V/4 



v 71 ' 



' 216 ik 3 / 2 ' 


1/2 - 


_sin 3 (VH)_ 





p- 



QiVk 
tan(Vkt)_ 



i -3/2 



( 6iVk 
x R exp{ 



1 + 



Qi^/k 



ltan(v / fct) V cos 2 {y r kt)[(3tan(^/kt) - 6iy/k}_ 



R 2 



(4.14) 



(4.15) 



The expectation value of the scale factor is found to be 

= 2 <£)i a) , (4.16) 



so that there is no singularity for boundary condition (3.9b) either. 

As a matter of fact, we have shown only that the states evolving from (4.1) or (4.13) are such 
that (R)t never vanishes. Incomplete analyses like ours of the quantum gravitational collapse 
problem have been made before [7,9]. This is insufficient, however, because in order to establish 
that the quantum cosmological models are nonsingular one has to prove that (R) t ^ for any 
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evolving state ip(t) for which (R) t is defined. A somewhat indirect proof of this will be given 
below. 

Classically the presence of the singularity at t — makes it physically mandatory to restrict 
the conformal time t to positive values. The absence of singularity at the quantum level makes 
such a restriction unnecessary, so that — oo < t < oo and the quantum cosmological models are not 
naturally endowed with an origin of time. This kind of situation is also encountered in dust-filled 
FRW models in the cosmic-time gauge [11]. 

4.2 Special Initial Condition 

Having a quantum dynamical framework to describe the evolution of the universe is not enough 
to explain its present state, one has to face the problem of initial conditions, the gist of mod- 
ern quantum cosmology. In the path-integral approach to quantum cosmology both the Hartle- 
Hawking and the Vilenkin-Linde proposals appear to suffer from vagueness and lack of generality, 
and can hardly be said to lead unambiguously to a unique universal wave function [29]. In the 
present minisuperspace model the Vilenkin-Linde tunneling boundary condition is not even imple- 
mentable, as we have seen in Section 3. 

With this in mind, we proceed tentatively to examine the outcome of imposing initial conditions 
on the probability density associated with the universal wave function. We content ourselves with 
discussing the extreme situation a — > oo, in which case \ipo\ 2 becomes sharply concentrated around 
R = : 

Jim \4 a) (R)\ 2 = S(R) . (4.17) 

Under such circumstances the universe starts with certainty from the singularity R = 0, so that 
(4.17) may be regarded as the condition for a quantum explosive birth of the universe, or what 
might be called a quantum big bang. For a sufficiently large Eq.(4.5) reduces to 

sint if k — +1 

t iffc = (4.18) 
sinht if k — — 1 

so that the classical regime sets in — compare the above equation with (2.10). 

Now for the promised proof that (R) t ^ for any evolving state ip(t), implying that our 
quantum cosmological models are nonsingular. If (R) tl = for some t 1 then \ip(R, ti)\ 2 = 5(R). 
Suppose that ip(R,t) with t > t± is an state evolved from ip(R,ti) taken as initial condition. By 
letting a — > oo it follows from (4.18) that (R) t = oo. Therefore, no state with finite expectation 
value of the scale factor can arise from ip(R,ti). Since quantum mechanics is time reversible, no 
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state ip(R,to) with t < h and finite {R)t can evolve to ip(R,ti), which proves that (R) t ^ for 
all evolving states ip(t) for which the expectation value of the scale factor is finite. 

Let us focus our attention particularly on the closed model (k=+l), the only one for which the 
following considerations are meaningful. In our present treatment, that deals only with normalized 
wave functions, the probability P(R < Ri,t) that at time t the radius of the universe is smaller 
than a given radius R\ is given by 

P(R<R i; t)^ = [ Rl \^ a \R,t)\ 2 dR = (—) 

JO V 7T ' 



36 

cos 2 1 [a 2 tan 2 1 + (6 — p tan t) s 



X io e ^\-cosn[aH^t + (6- P t^m} dR - (419) 



With the change of variable 



' = - — • - " ^ ^ 11/2 R , (4-20) 



cos t[a 2 tan 2 1 + (6 - ptant) 2 ] 1 ^ 



one gets 



P(R < Ri, t) (a) = ^= [ R * e~ x2 dx (4.21) 

\/7T JO 



where 



D = D (A 02) 

* cost[a 2 tan 2 t + (6-ptant) 2 ] 1 /2 /Tl • ^ > 



Notice that -R* — > as o" — > oo , hence P(i? < t) = for an explosive quantum beginning of 
the universe. This leads to the prediction that the density parameter Q equals unity in the limit of 
a truly explosive birth of the universe. Such a prediction was called "inflation without inflation" 
by Tipler [8], but here it is derived without having to resort to non-normalizable wave functions. 

If the initial wave function is chosen as (4.13), corresponding to the boundary condition (3.9b), 
Eq.(4.16 ) shows that if a is sufficiently large the classical regime takes over. For the closed model 
one readily finds 

P(R < Ri, t) (6) = -4= j R ' x 2 e~ x2 dx . (4.23) 
Vvr Jo 

Again in the limit a — > oo the probability that at time t the radius of the universe is smaller than 
a given radius R\ is zero. Note that as a — > oo the wave packet (4.13) also satisfies (4.17), thus 
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providing a concrete illustration of the fact that an initial wave function whose associated proba- 
bility density is concentrated entirely at R = can be harmonized with the boundary condition 
(3.9b). Besides, the prediction Q = 1 does not appear to be sensitive to boundary conditions of 
the type (3.6) on the wave function itself, but to result exclusively from the condition (4.17) that 
its modulus squared be sharply concentrated at R — when t = 0. An initial condition on the 
wave function itself that gives rise to Q = 1 is known [8], but then one has to give up the square 
integrability requirement, and the resulting universal wave function does not belong to the domain 
of self-adjointness of the Hamiltonian operator, as remarked in Section 3. 

Unfortunately, this state of affairs is still physically dubious. The inevitable singularity makes 
the restriction t > dynamically obligatory in classical cosmology. On the other hand, it is only 
the imposition of the initial condition \ip (R)\ 2 = S(R) that makes the instant t = so especially 
distinguished as to induce the restriction t > upon the quantized model too. This is so because 
no unitary evolution from an earlier time could have led to such a perfectly localized state at t — 0. 
This questionable feature is also present in Tipler's treatment [8], in which use is made of non- 
normalizable wave functions. Although inconclusive, our tentative considerations were intended 
to suggest that it may be physically reasonable to impose initial conditions on some probability 
distribution engendered by the wave function rather than on the universal wave function itself. 



5. QUANTUM STATIONARY GEOMETRIES 

In the case of the closed FRW model the Hamiltonian operator (3.1) possesses square-integrable 
eigenfunctions. Consider the normalized wave functions 

/ /To \ 1/2 

MR) = [ Jf 2 n n[ ) H - (SM R ) exp(-6i? 2 ) (5.1) 
where H n denotes the n-th Hermite polynomial. They satisfy 

H V>n = (™ + 1/2) (5.2) 



where n is a non-negative integer. For even n the wave functions (5.1) satisfy both Eq.(5.2) and the 
boundary condition (3.9a), whereas for odd n they obey both Eq.(5.2) and the boundary condition 
(3.9b). The effect of the boundary conditions is to exclude either even or odd eigenfunctions. 
This does not agree with [9], where no boundary conditions are imposed on the wave functions 
belonging to the domain of H, with the result that all values for n are allowed. For an universe 
in any of the stationary states (5.1) nothing changes with time. This is a purely quantum effect 
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since radiation-filled FRW universes do not possess classical static solutions. 

A more significant effect of imposing self-adjointness boundary conditions is the preclusion of 
oscillating coherent states, that is, nondispersive Gaussian wave packets whose center oscillates just 
like a solution to the classical equations of motion. From the general form of such wave functions 
[31] one recognizes at once that the general boundary condition (3.6) cannot be satisfied even if 
the real parameter a is allowed to be time dependent. It is clear, therefore, that a coherent state 
such that the classical solution (2.10) emerges as the expectation value of the scale factor does not 
exist, in contradiction with the findings in [9]. 

6. CONCLUDING REMARKS 

In this paper we dealt with square-integrable wave functions only, that is, we limited ourselves 
to the orthodox framework of quantum mechanics in Hilbert space. Interpretational controversies 
apart, it is in this arena that the requirement of self-adjointness on the quantum observables is most 
naturally justified and easily understood. Accordingly, being careful about domains of operators 
becomes a necessity in quantum cosmology, and in the case of radiation-filled FRW universes the 
simplest boundary conditions required to enforce self-adjointness of H were taken into account. 
Initial or boundary conditions introduced with the purpose of selecting a unique wave function are 
customarily unrelated to the former, although on occasion this has been object of confusion in the 
literature. Sometimes, however, these two types of boundary conditions interfere with each other, 
as we pointed out in Section 3. 

An important distinction should be emphasized between the classical and quantum cosmologies 
discussed in this paper. An origin of time comes into being dynamically in classical gravity due 
to the inevitable singularity, while no origin of time occurs naturally in quantum gravity, except 
if induced by a choice of initial conditions. On the other hand, admitting the hypothesis that a 
suitable initial condition exists, it is open to doubt whether it should be imposed on the wave 
function itself or on some probability distribution derived thereof. 

Our treatment differs from the one based upon consideration of conformal fluctuations about a 
given geometry. Apart from lack of application of boundary conditions to ensure self-adjointness, 
there is a physically more important difference. In [9] fluctuations are discussed about an specific 
FRW classical solution of Einstein's equations, whereas the ADM approach considers fluctuations 
that encompass all possible universes of the FRW kind. Although these standpoints may be 
classically undistinguishable, they are not necessarily equivalent in the quantum realm. This 
raises the question what the relation between these approaches is and which, if any, is appropriate 
from the physical point of view, a study we reserve for the future. 

We have circumscribed our analysis to radiation as cause of curvature, that is, to a matter 
content consisting of a perfect fluid with polytropic index 7 = 4/3. It can be shown [32] that in the 
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conformal-time gauge the form of the classical equation of motion of all Friedmann models with any 
7-ffuid as source can be reduced to that of a harmonic oscillator after a suitable change of variables. 
This suggests the possibility of extending the previous quantum treatment to cosmological models 
whose matter content is a perfect fluid with an arbitrary polytropic index. This is presently under 
investigation. 
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